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Abstract We apply the Basis Light-Front Quantization (BLFQ) approach to the Hamiltonian field 
theory of Quantum Electrodynamics (QED) in free space. We solve for the mass eigenstates corre- 
sponding to an electron interacting with a single photon in light-front gauge. Based on the resulting 
non-perturbative ground state light-front amplitude we evaluate the electron anomalous magnetic mo- 
ment. The numerical results from extrapolating to the infinite basis limit reproduce the perturbative 
Schwinger result with relative deviation less than 0.6%. We report significant improvements over pre- 
vious works including the development of analytic methods for evaluating the vertex matrix elements 
of QED. 
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1 Introduction 



Hamiltonian light-front theory has been proposed as a promising method for solving bound states 
problems in quantum field theory a long time ago, see, e.g., the review article [T] for its various advan- 
tages compared to the traditional instant-form theories. Recently the Basis Light-Front Quantization 
(BLFQ) approach has been developed as a nonperturbative approach to solve Hamiltonian light- 
front quantum field theory. Numerical efficiency is a key advantage of this approach. The basic idea 
of BLFQ is to represent the theory in an optimal basis which respects many symmetries of the the- 
ory and thus minimizes the dimensionality of the Hamiltonian for a fixed precision. Specifically, the 
BLFQ approach employs a plane wave basis in the light-front longitudinal direction and a 2D har- 
monic oscillator basis in the transverse directions. In previous work [3] this approach has been applied 
to evaluate the anomalous magnetic moment of electrons which are confined in an external trap with 
an extrapolation to the zero trap limit. In this work we extend and improve this approach in several 
aspects including the direct evaluation of a free electron system. This article is organized as follows: 
In Sec. [5] we discuss the key extensions and improvements made in this work over Ref [3J; in Sec. |3J 
we present the numerical results for the electron anomalous magnetic moment evaluated in different 
harmonic oscillator bases and compare to the perturbation theory result. Finally we conclude and give 
an outline for future works in Sec. 2] 
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2 Extensions and Improvements 

As in the previous work [3] we apply BLFQ to QED and solve for the mass eigenstates for single 
electrons in a truncated basis including only |e) and \e^/) sectors. In this paper we briefly report the 
major extensions and improvements made over Ref. [3]. A more systematic description for applying the 
approach to QED will follow in a separate paper [3] . Before introducing the details of these extensions 
and improvements, we list them here. 

1. We extend the application of BLFQ to a free space system by omitting the external trap. 

2. We replace numerical integrations previously used to evaluate matrix elements of interaction 
vertices in the QED Hamiltonian with newly-developed analytic methods. 

3. We correct the evaluation of the anomalous magnetic moment from the electron light- front am- 
plitude and a factor in the vertex matrix elements. These corrections go in opposing directions for the 
external trap results for the anomalous magnetic moment in Ref. [3j . The results of Ref. [3] with an 
external trap will be fully updated in a future publication. 

In the previous work [3] , the system (a physical electron consisting of a constituent bare electron 
and a bare photon) was configured with an external trap provided by a scalar harmonic oscillator 
(HO) potential acting on both the (constituent) electron and photon in the transverse directions. The 
frequency oj of the HO basis was chosen to coincide with the frequency of the HO trap potential. The 
motivation for introducing this potential is two-fold: a) it provides an ultraviolet cutoff for the QED 
interactions; b) it simplifies the Hamiltonian matrix since the basis is made of the eigenstates of the 
nonintcracting Hamiltonian. However this external potential may be difficult to implement experimen- 
tally due to its scalar nature and the fact that it confines both electrons and photons in the same trap. 
Therefore in order to compute quantities which are directly comparable with trap-free experimental 
data, it is desirable to set up the BLFQ approach in free space. To achieve this goal we implement the 
following procedures in the order we now list. 

1) We decouple the HO frequency u> chosen for the basis from that of the trap potential. One of 
the limitations in Ref. [3] is that the oj for the external trap potential was taken to be the same as 
that for the HO basis. However, in order to calculate observables in free space the oj for the trap 
potential needs to be taken to zero, while the oj for the HO basis has to remain finite. Note that, in 
principle, the basis is arbitrary as long as it is complete. Thus arbitrary w's for the HO basis should 
lead to identical physical observables. However, different choices of oj will exhibit different rates of 
convergence with increasing basis dimension to the "true" value, due to different effective infrared and 
ultraviolet cutoffs associated with each oj. Furthermore, the rates of convergence will also depend on 
the specific observable evaluated. We will illustrate this point for the anomalous magnetic moment in 
Sec. [3J Decoupling these two w's increases the number of nonzero off-diagonal matrix elements of the 
Hamiltonian matrix, but only mildly, because the kinetic energy term in the Hamiltonian, ~ p 2 , in the 
HO basis only couples states with adjacent principle quantum numbers leading to tridiagonal kinetic 
components in the Hamiltonian. 

2) We remove the transverse center-of-mass (cm) kinetic energy from the Hamiltonian. For the free 
space states only the relative transverse motion of the constituent particles is relevant for its invariant 
mass spectrum and light-front amplitudes. However, since BLFQ is set up in the single-particle coor- 
dinates, the kinetic energy of the transverse cm of all the constituent particles is naturally included 
in the sum over kinetic energies for all particles. Therefore the explicit removal of the transverse cm 
kinetic energy from the Hamiltonian is necessary. We achieve this by making the following replacement 

for the transverse momenta of all constituent particles n in a given Fock sector: p. L — > p\ = p. L — ■ 
It is easy to verify that p' t = and therefore the p\ carries only the relative transverse momentum. 

3) We remove all the excited configurations and keep only the lowest configuration in the single- 
electron sector, |e). With the cm kinetic energy removed from the Hamiltonian, the basis is overcom- 
plete for representing the internal motion among the constituent particles. This overcompleteness may 
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lead to multiple copies of the mass eigenstates with identical internal motion but different cm motions, 
which cause unnecessary numerical burden and complicate the mass spectrum. Therefore, it is desirable 
to remove these (spurious) copies from the spectrum. In this work we achieve this by simply removing 
all the excited configurations in the single-electron sector. Since these configurations produce states 
with excited states of cm motion, the removal of these configurations does not affect the completeness 
of the basis for representing the internal motion of constituent particles through higher Fock sectors. 
Furthermore, since the single-electron sector is the dominant sector for (physical) electrons in QED, 
removing the (spurious) copies of the cm motion from the single-electron sector effectively eliminates 
cm excitations from the lowest mass eigenstates. That is, the Hamiltonian of QED conserves cm motion 
and all excited states of cm motion will now involve only the \ej) configurations not coupled by QED 
to the |e) configuration that dominates the lowest mass eigenstate. 

4) We implement sector-dependent HO basis frequencies (w's). Different Fock sectors are coupled 
by the interaction terms in the QED Hamiltonian. These interaction terms are identity operators with 
respect to the transverse cm motion. Now, in the single-electron sector we already chose the cm mo- 
tion to be in the lowest allowed HO state. Thus, it is desirable to adopt the same state of cm motion 
for all higher Fock sectors. For Fock sectors with n particles, the cm motion can be described by 
HO bases with the characteristic length parameter b n = y/nMuj n , where M is the mass of the HO 
specifying the basis. (Note that this M has nothing to do with the physical mass of particles in the 
QED Hamiltonian and, for convenience, all single-particle basis states in all Fock sectors are assumed 
to have the same M). In order to match b n for all Fock sectors, the HO frequency in each n-particle 
Fock sector is required to satisfy 0J n = h UJ i) where u>i is the basis frequency for the single electron sector. 

The second improvement is related to the calculation of the matrix elements of the interaction 
vertices between different BLFQ states. Let us take the matrix elements coupling |e) and \ej) sectors 
for example. To evaluate these matrix elements one needs to perform integrations in the transverse 
momentum space over a product of three HO wave functions associated with the e in |e) and the 
e,7 in |cy} sectors. In Ref. [3] these integrations were computed numerically in the single-particle 
representation. Due to the oscillatory structure of HO wave functions, the numerical computation was 
very time-consuming to achieve results with the desired precision. In this work wc reduce the numerical 
integrations by applying an analytical Talmi-Moshinsky transform [7] to simplify the integrations. The 
needed 2D Talmi-Moshinsky coefficients are analytically worked out in Ref. [8]. 

Adopting this analytic step greatly reduces the calculation time for the desired precision of the 
resulting matrix elements. For example, numerically integrating all matrix elements coupling |e) and 
\ej) sectors with the truncation at iV max =4C0 previously took ^7000 cpu-hours at the National Energy 
Research Supercomputer Center (NERSC) and only ^15% of them have relative errors better than 
one part in I0 7 . With the improved methods, we now evaluate these same matrix elements in less 
than 2 cpu-seconds and all of them have relative errors better than one part in 10 7 . Developing and 
adopting these analytic techniques for interaction matrix elements solves a bottleneck in BLFQ and 
greatly expands the attainable basis space dimensions. Increasing basis space dimension is crucial for 
evaluating observables with slow convergence with N max and/or improving the overall precision. We 
expect that the same analytic techniques can be applied to calculate all interacting vertices in QED 
and QCD. The details for the analytic calculation of the interaction matrix elements will be given in 
a forthcoming paper [¥]. 

Our third improvement is in the method of extracting the electron anomalous magnetic moment 
from its light-front amplitude obtained from diagonalization of the QED Hamiltonian matrix. We also 
correct a factor in the calculation of the vertex matrix elements in the Hamiltonian matrix. These two 
corrections to the results in Ref. [3] have approximately offsetting consequences. 

In Ref. [3] the electron anomalous magnetic moment was approximated (based on non-relativistic 
quantum mechanics) by the squared modulus of the helicity-flip (for the constituent electron) compo- 
nents of the eigenstates. The precise definition of electron anomalous magnetic moment in relativistic 



1 In BLFQ TVmax indicates the upper limit for total HO quantum numbers summed over all particles in a 
Fock sector. Any state whose total HO quantum number exceeds JV mm are omitted. For -/V max =40 there are 
8190 unique non-vanishing electron to electron-photon vertex matrix elements in total. 
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QED is a e , the electron Pauli form factor F 2 evaluated at momentum transfer q 2 — > [6], 

a e ^^=F 2 (q 2 ^0). (1) 

In BLFQ the a e can be calculated by sandwiching the operator ^2(0) with the obtained ground states 
for the electron, 

ae = ^ = <^|F a (0)|ff-t) 

= J2{^\e 1 ,i'}(e 1 ,i'\F 2 (0)\e 1 ,i){e 1 ,i\^). (2) 

i' ,i 

Here (ey, i'l-F^O)^, i) is the matrix element of the Pauli form factor in the BLFQ basis. The (ej, ) 
is the wavefunction of a physical electron with helicity up (down) in the \ej) sector (the only sector 
contributing to a e in our truncated basis). The i denotes a complete set of quantum numbers. Although 
Eq.@ involves two electron eigenstates with opposite helicities, in practice one needs only to solve 
for one of them and infer the other by exploiting the parity symmetry in light-front QED [51. The 
explicit expression for (ej, i'\F2(0)\ej, i) and the exact relation between (cy, zltf^J") and (e~f, will 
be reported in a later work [4|. 

As noted above, we also correct the interaction matrix elements used in Ref. [3] where there was an 
inadvertent omission of a factor of 2. This, in turn led to an overestimate of the evaluated square root of 
anomalous magnetic moment based on the spin- flip modulus by a factor of 2. Fortuitously, according to 
an analysis based on the perturbation theory 0] , the modulus of the helicity flip amplitude is down by 
a factor of y2 from a e defined above. Therefore, these two corrections are partially offsetting. Further 
improvements to the methods of Ref. [5] for the extrapolation to the free space limit are now envisioned 
based on the sensitivities we observe in the results presented here. The revised results for an electron 
in an external trap will be presented in a later work [3] and will include results in significantly larger 
basis spaces to underpin an improved extrapolation to the free-space limit. 



3 Numerical Results for Electron a e 

Incorporating all the improvements introduced in Sec. [5J we diagonalize the light-front QED Hamilto- 
nian in |e) and jey) sectors and extract the anomalous magnetic moment from the resulting light- front 
amplitude for the lowest mass eigenstate. In this work we reduce the QED coupling constant a by a 
factor of 10 7 to limit the difference between our non-perturbative result and the perturbative result. 
Then, in this small coupling limit, the contributions from higher order perturbations can be neglected 
so that the Schwinger result [5] can be used as the reference for our numerical results. 

We evaluate and diagonalize the Hamiltonian with total longitudinal momentum (summed over the 
constituent electron and photon) _ftT=160, translating to 160 (equal distance) partitions between the 
longitudinal momentum of the constituent electron and photon. In the transverse directions we use HO 
basis with several selected frequencies cj's, ranging from O.OlMeV to 1.4MeV bracketing the electron 
mass m e =0.511MeV, the only dimensionful quantity in the QED Hamiltonian. At each oj we calculate 
a series of results with increasing N max from 10 to 118, for the purpose of mapping out the convergence 
behavior of anomalous magnetic moment with respect to N max . The larger iV max translates to a larger 
basis (the largest here consists of 44,358,561 basis states) with higher effective ultraviolet cutoff and 
lower effective infrared cutoff in the transverse plane. It is thus expected that, with increasing iV max , 
the results grow closer to the Schwinger result. The rate of convergence may, however, be different 
for different w's, depending on the sensitivity of the anomalous magnetic moment to the infrared and 
ultraviolet cutoffs in the transverse directions. Indeed our numerical results agree with this expectation: 
all the results with w between O.OlMeV to 1.4MeV are below the Schwinger result at finite -/V max which 
they approach uniformly as -/V max increases with increments of 4. 

In Fig. [T] we compare the results evaluated with w=0.02MeV and w=0.5MeV. At the same N max 
the results at w=0.5MeV are closer to the Schwinger result than those at w=0.02MeV. For each 
oj the results exhibit a distinctive pattern with increasing -/V max : the results with even N max /2 are 
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Fig. 1 (Color online) Anomalous magnetic moment of the electron calculated in BLFQ compared to the 
Schwinger result [9] . The vertical axis is the square root of anomalous magnetic moment normalized to electron 

charge, e, so the Schwinger value is \J~^i = 0.11254. The horizontal axis is the square root of the reciprocal 

of iV max . Symbols are for the BLFQ results. Squares: cj=0.5MeV at even 7V max /2; circles: ct>=0.5MeV at odd 
Af m ax/2; triangles: aj=0.02MeV at even -/V max /2; diamonds: ct>=0.02MeV at odd iV max /2. The lines are linear 
extrapolations of BLFQ results based on the points at or above iV max = 64. 



systematically larger than those with odd N max /2 so that the former and the latter separate into two 
individual groups. Within each group the results define a simple trend. 

The systematics shown in Fig. Q] have simple explanations. According to perturbation theory [B], 
the a e receives its largest contributions from the infrared region of the relative momentum between the 
constituent e and 7. Thus, a difference between neighboring even and odd N max points originates from 
a truncation-induced distortion of the (point-like) contact interactions (coupling |e) and \ej) sectors) 
in the infrared region of momentum space. If one works in the relative coordinates (for the constituent 
e and 7), representing such (point-like) contact potentials in a truncated HO basis typically leads 
to an alternating over- and under-estimate in the infrared region of the relative momentum at even 
(odd) N mayi /2 truncations. It is natural to expect that such qualitative features are the same in the 
single-particle coordinates used in BLFQ. 

The uj dependence of the convergence is also understandable. Since the HO basis with smaller 
u/s have more accurate representation for the infrared region, the odd-even oscillation of results at 
w=0.02MeV is milder than those at w=0.5MeV. On the other hand, the results at w=0.5MeV are 
closer to convergence at any chosen N max reflecting naturalness of this basis for the QED problem 
with only the electron mass scale. 

For sufficiently large N meiK the data points in Fig. [1] define straight lines as a function of l/y/N max . 
We can therefore easily extrapolate to the limit of without basis truncation (N max — > 00) where we 
expect to recover the Schwinger result. Indeed as seen in Fig. Q] the four extrapolating lines con- 
verge close to the Schwinger result in this limit. Their intercepts at 1/N max —Q are: 0.11285(0.5%), 
0.11248(0.1%), 0.11255(0.02%) and 0.11280(0.5%) for w=0.5MeV at even N max /2, w=0.5MeV at 
odd N max /2, w=0.02MeV at even N max /2, w=0.02MeV at odd N max /2, respectively. The percent- 
ages in the parenthesis are their corresponding relative deviation from the Schwinger result, % = 
7T^z=^K- »0. 012665. These results suggest that the u/s between 0.02MeV and 0.5MeV form an "op- 
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timal" plateau, where the Schwinger result is reproduced with less than 0.6% deviation. We also note 
the large width of this plateau (from u;=0.5MeV to 0.02MeV amounting to a factor of 25), illustrating 
the flexibility and stability of the BLFQ approach. For u> slightly below 0.02MeV or above 0.5MeV 
reasonable precision can still be achieved and details will be presented in a later work 0]. Finally, we 
would like to emphasize the importance of large V max results in extracting high precision a e , since 
only at sufficiently large values of N max do these points follow a straight line in 1/ V-^max, allowing for 
reliable linear extrapolations as seen in Fig. [T] 

4 Conclusions and Outlook 

Following Ref. [3] we apply the BLFQ approach to QED and solve for the electron anomalous magnetic 
moment. Several extensions and improvements are introduced over the previous work By removing 
the transverse external trap we are now able to solve for the mass eigenstates of QED in free space 
so that the resulting observables are directly comparable with most experimental data. The analytic 
procedures introduced for interaction matrix elements greatly reduce the total computing time and 
improve the overall precision. As a result, we can now access much larger basis spaces in BLFQ calcula- 
tions, which is crucial for evaluating observables having slower convergence rates. As an application of 
BLFQ we evaluated the electron anomalous magnetic moment and compared to the Schwinger result. 
We found that the extrapolated values agree with the Schwinger result within reasonable precision 
with the present truncations. We found that high precision results can be obtained with HO basis 
spaces over a wide range of frequencies, signifying the flexibility and stability of the BLFQ approach. 

We plan further developments. 1) We will investigate renormalization schemes for light-front field 
theory such as those presented in Refs. [TUlfn] and calculate QED with the physical coupling. 2) We 
will extend the basis space to include higher Fock sectors and include all the interacting vertices in 
the QED Hamiltonian (including instantaneous vertices). These steps will allow BLFQ to capture all 
the main physics in QED and make it an ideal approach for investigating QED in nonpcrturbative 
regimes, such as in the presence of strong pulsed laser fields. 3) Furthermore, we will apply the BLFQ 
approach to QCD, where anticipate the BLFQ approach will be able to provide valuable information 
on the nonperturbative aspects of QCD. That is, in one perspective, we will extend the AdS/QCD 
approach of Ref. [H] to include QCD interactions. 
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